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and density dependent particle mass is studied. Many previous treatments have been reviewed. 
A new thermodynamical treatment based on the equilibrium state is suggested. Employing the 
quark mass density- and temperature-dependent model, the discrepancies between our treatment 
and others are addressed. 
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I. INTRODUCTION 



It is generally accepted that the effective masses of particles will change with temperature 
and density due to medium effects. Many theoretical considerations, including the finite 
temperature QHD modeljlj], Brown- Rho scalingj^, QCD sum rules j^], finite temperature 
QMC model [4], |5j, vacuum polarization Feynman diagrams calculations , 0] and etc., 
had been suggested to investigate the effective masses of mesons and nucleons. By means of 
the Thermo-field Dynamics, after summing the temperature- and density-dependent vacuum 
polarization diagrams for tt, cr, lu, p mesons and the three-line vertex diagrams which give the 
corrections of the NN71", NNcr, NNu^NNp couplings respectively, we have extended the one- 



9]. ()ii< , -])H)u|10] ! one-omega and one-rho exchange potentials 0, [g, Q, and then 
13| to finite temperature and density We have found that the attractive 



sigma 



the Bonn potentials 



part of the NN interaction becomes weaker and the repulsive part becomes stronger when 
the temperature and/or density of the system increase. This is of course very reasonable. 

Besides theoretical study, many experimental results which predict the changes of particle 
masses with temperature and density have been shown. In particular, the experiments of 
TAGX collaboration have shown directly that when the density of the nucleon medium 
equals to 0.7no where no is the saturation density, the effective mass of nutral p-meson 
reduces to 610Mev[l^]. Both theoretical and experimental results confirm that the medium 
effects are important for studying the nuclear or quark systems. 

To illustrate the medium effects more transparently in the theory instead of the first 
principle calculation, many authors introduced different hypothesis to represent the medium 
contributions to particle masses, for example, introduced the density-dependent vacuum 
energy B(n) to modify QMC model 15, 



to address liquid-gas phase transition 



16| . suggested the density-dependent NNp coupling 



171 ]. supposed the u, d and s quark masses depends on 



density to mimic the quark confinment mechanism [181]. introduced temperature-dependent 
vacuum energy B(T) to modify the quark mass density-dependent (QMDD) model 191 ]. and 
etc.. Employing these hypothesis, many physical properties of nuclear matter, quark matter, 
nucleon system and hyperon system had been discussed. 

Although the density- and temperature-dependent particle masses m*(T,n) can mimic 
the medium effects, when we discuss the thermodynamical behaviors of the system with 
such particles, many difficulties will emerge. First, the dispersion relation for a particle with 
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energy e and momentum k becomes 

e(k,T,n) = [k 2 + m*(T,n)} 1/2 (1) 



due to the medium effect. It was shown in Ref.j20| that the usual fundamental thermody- 
namical partial derivative relation is not satisfied for such system, we must add the additional 
terms. Secondly, it is well known in thermodynamics that if we choose a pair of independent 
variables, for example, temperature T and volume V, and give the suitable characteris- 
tic thermodynamical function, for example, Helmholtz free energy F = F(T,V), then all 
thermodynamical quantities can be obtained by the partial derivatives of F(T, V), because 

dF = -SdT - pdV, (2) 
OF dF 

S = -(gf)v, p=-{gy) T . (3) 

Eq.(3) give the entropy and the equation of state respectively. Other thermodynamical 
quantities such as internal energy U, enthalpy H, Gibbs function G, heat capacities cy and 
c p can be obtained from Eq.(2) by derivatives. In above calculations, no integral process 
will be needed and no integral constant has to be determined. In fact, this is the reason 
why all calculation of thermodynamical quantities in the canonical ensembles can be done 
with the calculation of the partition function Z = Z(T,V), because the partition function 
Z is related to free energy F directly. For systems with three independent variables, say, 
T, V and chemical potential /i, the suitable characteristic function is the thermodynamical 
potential Q = Q(T, V, n) because 

dVL = -SdT - pdV - Ndfi. (4) 

But for system with density-dependent mass particle, Q is not only a function of T, V, /i, but 
also of the density n, because Q = Q(T,V, /j,,m*(T,n)). How to treat the thermodynamics 
self-consistently is still a serious problem and has many wrangles in present references. We 
will give a brief review of various treatments and show their contradictory in the next section. 
A few comments on different treatments will also be presented there. 

This paper evolves from an attempt to suggest a new thermodynamical treatment to 
study the system with density- and temperature-dependent mass particles. Since the particle 
mass depends on temperature and density, a lot of ambiguities will happen when one uses 
the partial derivatives along a reversible process to obtain the thermodynamical quantities. 



To avoid this difficulty, instead of a reversible process, we argue that we can calculate the 
thermodynamical quantities at an equilibrium state. In fact, all the physical quantities, 
such as p, U, F, H, G, are the functions of equilibrium state and have definite values 
respectively. We can use the relations between thermodynamical functions to obtain above 
quantities provided that the thermodynamic potential Q can be obtained. We will employ 
the QMDD model to explain our treatment in detail in Sec. III. By means of quark mass 
density- and temperature-dependent (QMDTD) model, we will show our results and the 
comparison to other treatments in Sec. IV. The last section is a summary. 



II. DIFFERENT THERMODYNAMICAL TREATMENTS 



There have been several different kinds of treatments dealing with the thermodynamics for 
the system with temperature- and/or density- dependent mass particles. The basic difficulty 
comes from Eq.(l) while one calculate the partial derivative of thermodynamical functions 
along a reversible process. Though several authors had still employed the usual direct 



formulae to calculate the thermodynamical quantities and shown 



21 



22 



P = - V , (5) 

i 

where ub is the baryon density, rii is the number density of particle i, \ii is the corresponding 
chemical potential, and e is the energy density. Their failure is obviously because, with the 
temperature- and density-dependent mass term, Q is no longer an explicit function of T, 
V, Hi. The usual direct partial derivative formulae are not applicable. To overcome this 
difficulty, many different treatments or methods have been suggested in the market, which 
can be generally categorized into two kinds. 

(I) The first kind of treatments takes the mass term as a function of density and/or tem- 
perature and make a great effort to work through the derivatives of a composite function 
along a reversible process. In deriving the energy density and the pressure, extra terms 
emerge in this approach due to the dependence of the mass on the density and/or temper- 
ature. Although many authors work along this direction to study the mass density and/or 
temperature dependence, their formulae and results are very different with each other. 



(IA). For the the density- dependent mass, in Ref. [231] . they gave 



1 d{tt/n B ) tt n B dtt 

n n B .dn. t m. 

where rij still satisfies Eq.(6). The extra terms in Eqs.(8, 9) produce significant changes in 
the energy density e, and makes that the pressure become negative in low density regions. 



(IB). In Ref. 24], |25|, for the density-dependent mass, the pressure and the energy density 
become 

_. tt n B , dtt , . 

p = -v + vtt™> (10) 
q t on 

S = y+l^^ n i- y(-Qf)UH},n B - (H) 

i 

They use the same pressure formula as that of (IA) Eq.(8), but do not agree with the 
expression in Eq.(9), because it cannot give a correct QCD vacuum energy. Then they use 
the same energy density formula as Eq.(7). 



(IC) . In Ref. |26j , for the density- and temperature-dependent mass, the pressure and the 



energy density read 

„ ~ T r dtt ^ dtt drrii , . 

p = - a ' v av +nB ^^ (12) 

~ ^ dtt dtt dtt drrii 

£=n -^^r T df' T ^d^^ (13) 

where tt = tt/V stands for the density of the thermodynamical potential. 

(II) The second kind of treatments add different terms to the system to keep the usual 
thermodynamical derivative relations unchanged. For example, in order to make the Eq. (4) 
and the calculation of the thermodynamical potential self-consistently, they add extra terms 
to the Hamiltonian or to the thermodynamical potential. These additional terms differ for 
different authors. 

(IIA). To keep the fundamental thermodynamic relations 

self-consistently, in Ref. a term Eq had been added to the Hamiltonian. This term is 
determined by the condition 

dp 

{g-)T^{ Cj ^} = 0, (15) 



where {q} is the temperature- and density- dependent terms in the Hamiltonian. Then the 
pressure and energy density become 

_gkT 



p(T, ^{ Ci }) = T 9 —j k 2 dk ln[l T exp-^*^] - B* 

e{k)k 2 dk 



e(T,fi,{ci}) 



where 
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2vr 2 J exp-P« k )-ri =pl 



B* = lim % 



B\ 



(16) 
(17) 

(18) 



and g is the degeneracy factor. 

(IIB). In Ref.j^, an extra term fi (ng) had been added to the thermodynamical poten- 
tial, 



9iT 
(2tt) 3 



dfc 3 ln(l + e ~/W))-w) + n a (n B ), 



where fi a (n B ) is determined by the constraint 



0. 



At zero temperature, the corresponding thermodynamical formulae become 

P = -0, 



(19) 



(20) 



(21) 



£ = Et^^W -m?(nfl)] 1/2 [6/i?-3m?(nB)] 



48tt 



-3m,> B ) ln[ "- + ^ - "fl"")! 1 * ]} + !!.(»,), (22) 



67T' 



rrii(n B ) 
•\A ~ m ii n B)} 3/2 - 



(23) 



In summary, from above treatments we come to a conclusion that how to treat the 
thermodynamics with the medium effect is still a serious problem and has made many 



wrangles in present references 
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23, 



24 



25|, 



26J. In fact, many treatments are contradict 



with each other. To show their confusion and ambiguity, we give a few comments in the 
following: 

(1). The first argument is that the dispersion relation Eq.(l) and the derivative for- 
mula of thermodynamical potential Eq.(4) correspond to different conditions. In Eq.(l), 
the medium effect has been taken into account. This effect makes that the particle mass 



becomes a temperature- and density-dependent effective mass. But Eq.(4) corresponds to 
the thermodynamical system only where fl, S, T, p are the thermodynamical quantities of 
the system respectively. They do not include the medium effect. In order to make Eq.(l) 
and the derivative formula of thermodynamical process self-consistent, instead of Eq.(4), 
we must establish a formula of thermodynamical process for the total system. It includes 
the variations of the quantities not only for the system, but also for the medium. In this 
formula, all extensive quantities must be replaced by the quantity of the total system, for 
example, S — > S s + S m , where S s and S m are the entropy of the system and the medium 
respectively. The problem is that so far we do not know how to calculate S m and other 
extensive thermodynamical quantities of the medium directly. 

(2) . To illustrate above argument more transparently, we imagine two systems. These 
two systems are almost identical: their only difference is that one system consists of constant 
mass (mi) particles while the other of temperature- and density-dependent mass (m 2 (T, n)) 
particles. Fixed the temperature and density as T and no respectively, which satisfy 

m 2 (T ,n ) =mi, (24) 

at the equilibrium state with (T ,no), obviously, the two system are completely the same. 
They arrive at and stay in the same equilibrium state and have the same thermodynamical 
quantities, such as the pressure and the energy density, which are just functions of equilib- 
rium state. But if we use equations with partial derivatives to calculate these quantities of 
the later system, extra terms, such as J^:§^ and J^-^ in Eqs.(12, 13), emerge. These 
terms will never appear in the calculation of the former system since in this system mi is a 
constant. This simple example tells us that the first kind of treatments is not right because 
the medium effect has been neglected in Eq.(4). 

(3) . Now we hope to give a brief comment on the second kind of treatments. They 
hope to add an extra term to the system to consider the medium effect. But unfortunately, 
in general, the extra term cannot be determined by an additional constraint in terms of 
partial derivatives uniquely. It cannot be expressed as a zero value of thermodynamical 
functions, because the effective mass m* (T, n) comes from the interaction of the particle 
and the medium. This interaction depends on temperature and density. As an example, 
let's discuss the treatment (IIB). Obviously, if we add an arbitrary temperature function 
f(T) to ila(ns), the new expression for total still satisfies the constraint Eq.(20), but the 
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additional arbitrary function of T will change the thermodynamical functions, such as the 
entropy, which depend on the derivation of temperature. 



III. NEW TREATMENT BASED ON THE EQUILIBRIUM STATE 

The above arguments impress us to give up the derivative calculation of thermodynamical 
quantities along a reversible process, since Eq.(4) dose not include the medium effect. Instead 
of studying the reversible process, we focus our attention on equilibrium states. We suggest 
a new treatment which is based on the equilibrium state in this section. We will show there 
is no ambiguity in our treatment. 

According to the thermodynamics of equilibrium state, all thermodynamical functions 
such as Q, S, U, G, p ... have definite values respectively at an equilibrium state. If we can 
calculate the thermodynamical potential Q and the Gibbs' function G for a fixed equilibrium 
state, we can use the following definitions to find other thermodynamical quantities: 

U = £n iei , (25) 

i 

n = — = y^n,- = — -, — , (+ for Fermion and — for Boson) (26) 

G = (27) 

i 

F = Q + G, (28) 
S = (29, 

P = (30) 



To show our treatment explicitly, we employ the QMDTD model [13, |28j, |29|, l30l l3ll] . The 
QMDTD model is extended from the QMDD model which was first suggested by Fowler et 
al. many years ago [18]. According to the QMDD model, the masses of u, d quarks and s 
quark are given by 

m g = - — , (q = u,u,d,d), (31) 
6ns 

B . 
m s ,s = m s0 + - — , (32) 
3n B 

where B is the vacuum energy density and m S Q is the current mass of the strange quark. It is 
clear that the masses of quarks become infinity when the baryon density goes to zero, which 
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means that the quark confinement in this model is permanent. To remove the permanent 



confinement of quark, we modified the QMDTD model by introducing 

B{T) 



m 



3n 



B 



m s0 + 



(q = u, u, d, d), 
B{T) 



where 



3n B ' 

B(T) = B [l - (^) 2 ] 



(33) 
(34) 

(35) 



The quark mass depends on both density and temperature, so it is called a QMDTD model. 



We have used this model to study the properties of strange quark matter (SQM) 



28 



30 



3l|. 



For the system of SQM, in which the masses of quarks satis 



state, the thermodynamical potential of the system reads 



23 



E 



QjTV 
(2tt) 3 



d 3 k ln(e 



qs. (33-35), at equilibrium 
26|,|28j 

(36) 



+ 1) 



The internal energy U and the Gibbs' function G are 



u -T,9iJ<Pk- e _ /3iei{k) _ iH) + l . 



(37) 
(38) 



e -p(ei(k)-iH) -f 1' 

respectively. Other thermodynamical quantities can be obtained by Eqs. (25-30). 

At finite temperature, the antiquarks must be considered. The baryon density satisfies 



where 

An,; 



rii - nj 



(27T) 



n B = -(An w + An d + An s ), 



d 3 k( 



exp[f3(ei(k) - fa)] + 1 exp[/3(e;(/c) + fa)} + V 



nj)rii is the number density of the (anti)flavor i{i = u, s, d), gi — 6, for antiquark fij 



Following Ref. [22 



28|, the system of SQM must satisfy the constraints 

Us = Hd = V-u + f*>e, 

because inside SQM, s (and s) quarks are produced through the weak process 

u + d u + s, s — > u + e~ + I7 e , d — > u + e _ + I7 e , u + e~ — * d + u e , 
and similarly for antiquarks. The condition of charge neutrality reads 

2An u = Ana + An s + 3An e . 



(39) 



(40) 



(41) 



(42) 



(43) 
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IV. RESULTS AND DISCUSSION 



For the convenience of comparison, our numerical calculations have been done by adoptin 



the parameters B = 170MeVfm" 3 , m s0 = 150MeV and T c = 170MeV, as that of Ref. 



28] 



Our results are summarized in Figs. 1-4 and Table I. 

The energy per baryon vs. baryon number density is shown in Fig.l where the 
temperature is fixed at 50MeV. In Fig.l, the solid line refers to the present treatment and 
other four dashed lines refer to treatments in Ref. 



21 



22j, IA, IB and IC respectively. We 



see the solid line is lower than the others. The saturation points for different treatments are 
summarized in Table I. We find from Fig.l and Table I that the differences are remarkable 
for different treatments. 

To illustrate the thermodynamical characters of our treatment, we show the internal 
energy and the free energy curves calculated by our treatment in Fig. 2 respectively. The 
energy per baryon ejn^ (solid lines) and the free energy per baryon //% (dashed lines) vs. 
baryon density for different temperatures T = 0, 30, 50, 100 MeV are shown in Fig. 2. We 
find that e/% increases and //% decreases as the temperature increases. The saturation 
points exist for all temperatures. 

The equation of state for different treatments are shown in Fig. 3, where the solid line 
refers to present treatment and other lines for previous treatments, respectively, as labeled 
in the figure. The solid line exhibits a significant property which differs from that of the lines 
of treatments IA, IB and IC. The pressure is definitely positive in our treatment. But for IA, 
IB and IC treatments, the pressure becomes negative in the small energy density regions. 
The pressure will not be negative due to its thermodynamical treatment for a system with 
positive energy. This result confirms that our treatment is correct. 



treatments 


riBo(fm~ 3 ) 


e/n B 


in refs. [21,22] 


0.46 


1023 


IA 


0.55 


1083 


IB 


0.36 


1120 


IC 


0.34 


1084 


our new treatment 


0.45 


1007 



TABLE I: Saturation points at T = 50MeV for different treatments. 
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In Ref. 



26|, the authors claimed that the dispersion relation of effective mass m* must 



satisfy the constraint 

Otti* 

lim — — = 0, (44) 



based on their formula for entropy 



dVL ^ dVL drrii 



If linij^o ^pf~ 7^ 0, it will conflict with the third law of thermodynamics, lim^^ 5 = 0. We 
hope to point out that this argument is not right because it is based on Eq. (4) which has not 
taken the medium effect into account. In Fig.4 we draw the entropy per baryon vs. temper- 
ature curves for different baryon densities n B = 0.2,0.5, 1.0 fm -3 with dispersion relations 
77i* (T, rie) = ^3^"^ ) using our treatment. We even choose three different dispersion re- 
lations m*(T,n B ) = Bo(1 3^ /Tc) , m* 2 (T,n B ) = Bo{1 ^ /T * ] and m* 3 (T,n B ) = £^ for the same 
baryon density n B = 0.5 fm~ 3 and find that all curves get together at the point T = 0, 
5 = 0. It means that our treatment is consistent with the third law of thermodynamics no 
matter how the dispersion relation is. The constraint Eq.(44) needs not to be satisfied. 



V. SUMMARY 



In summary, we have shown the shortcomings of the previous treatments, which based 
on the partial derivatives of thermodynamical functions along a reversible process or based 
on the additional terms to thermodynamical potential. The previous treatments obstruct 
the correct consideration of medium in thermodynamics. A new treatment of medium 
thermodynamics based on equilibrium state is suggested. Employing QMDTD model, we 
address the discrepancies between our treatment and others in Fig. 1-4 and Table I. We 
find that the negative pressure and the constraint condition Eq.(44) are removed in our 
treatment. 
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FIG. 1: Energy per baryon as a function of the baryon density ub at T = 50 MeV for different 
treatments [281 ]. The treatment IC and our treatment are for the QMDTD model, and others are 
for the QMDD model. 



14 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 



n B (fm" ) 

FIG. 2: Energy per baryon U/Nb = s/ns and free energy per baryon F/Nb = f /ns as functions 
of the baryon density ns at different temperatures T = 0,30,50,100 MeV. The two curves are 
identical at T = 0. 
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FIG. 3: Pressure P as a function of the energy density U/V for different treatments of both 



the QMDD and the QMDTD models 2a]. The tendencies of the curves are similar at large 
energy density region, but at small energy density region, different treatments have quite different 
behaviors. In our treatment, the pressure never goes to negative. 
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FIG. 4: Entropy per baryon vs. temperature curves for densities ns = 0.2,0.5,1.0 fm -3 , where 
the mass has been taken as m\{T, ng), m^T, ng) and m^T, ng) respectively, as indicated in the 
figure. All curves get together at S = when T = 0, which is consistent with the third law of 
thermodynamics. 
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